Theory of Thermoelectric Power in High-Tc Superconductors 
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We present a microscopic theory for the thermoelectric power (TEP) in high-Tc cuprates. Based 
on the general expression for the TEP, we perform the calculation of the TEP for a square lattice 
Hubbard model including all the vertex corrections necessary to satisfy the conservation laws. In 
the present study, characteristic anomalous temperature and doping dependences of the TEP in 
high-Tc cuprates, which have been a long-standing problem of high-Tc cuprates, are well reproduced 
for both hole- and electron-doped systems, except for the heavily under-doped case. According to 
the present analysis, the strong momentum and energy dependences of the self-energy due to the 
strong antiferromagnetic fluctuations play an essential role in reproducing experimental anomalies 
of the TEP. 



The thermotransport phenomena in high-Tc cuprates 
exhibit various anomalous behaviors. In particular, the 
thermoelectric power (TEP), expressed by the Seebeck 
coefficient S, is known to show characteristic temperature 
and doping dependences above Tc, namely (i) dS/dT < 
in hole-doped compounds such as YBa2Cu307_2; 
(YBCO) or La2-^Sr^Cu04 (LSCO), and S is posi- 
tive in the under-doped systems at room temperatures 
|Q,||. (ii) S is negative in electron-doped compounds 
like Nd2-^Ce^Cu04 (NCCO) and LaPri_a;Ce:,Cu04 
(PCCO), and dS/dT > in the under-doped compound 
at higher temperatures T > 200 K |g,|[. In both cases, 
|5| increases drastically as the doping decreases. Thus, 
a conventional Fermi liquid type behavior, 5 oc T, is 
totally violated in high-Tc cuprates for a wide range of 
temperatures. In particular, a qualitative particle-hole 
symmetric behavior of the TEP, i.e., ^lsco ~ — •S'ncco 
for the same carrier doping weight x as was reported in 
ref. ||], is very mysterious because both LSCO (YBCO) 
and NCCO have similar large hole-like Fermi surfaces 
according to angle-resolved photocmission spectroscopic 
(ARPES) studies |,|. 

From an academic point of view, the TEP is a unique 
and important phenomenon in that it sensitively re- 
flects the properties of the quasiparticles away from the 
Fermi surface, whereas electronic transport phenomena 
are caused by the quasiparticles on the Fermi surface. 
Thus, the characteristic non-Fermi liquid behavior of 
the TEP in high-Tc cuprates suggests that the excited 
states of quasiparticles in high-Tc cuprates are signifi- 
cantly anomalous. In this respect, the theoretical study 
of the TEP is quite important. 

Up to now, several theoretical studies on the TEP 
for models with strong Coulomb interactions have been 
performed based on the dynamical mean field theory 
(DMFT) [IP. These studies found that the energy 
dependence of the relaxation time is important for the 
TEP in strongly correlated systems. However, in high- 
Tc cuprates, it is known that the momentum dependence 
of the relaxation time is also prominent, and the vertex 
correction (VC) in terms of the current plays an essential 
role for transport phenomena. These features are totally 



dropped in DMFT. Thus, for the present purpose, we 
have to develop a theory for the TEP by taking the mo- 
mentum dependence into account. 

In this letter, we study the TEP in high-Tc cuprates 
on the basis of the antiferromagnetic (AF) fluctuation 
theory. We take the momentum and energy dependences 
of the relaxation time into account, and include all of 
the VC's required by the Ward identity, i.e., the con- 
serving approximation S. Our study reproduces the 
main futures of the experimental TEP for both hole- 
and electron-doped compounds. According to the present 
analysis, the approximate particle- hole symmetry for the 
TEP is caused by the " alternation of cold spots" , as is 
the case with the Hall coefficient [p|,[lO| 

In the present study, we calculate the self-energy Sk(e) 
using the fluctuation-exchange (FLEX) approximation 
which is a kind of self-consistent spin-fluctuation the- 
ory |11|], because it works well for high-Tc cuprates for 
|1 ~"l ^ 0-1 above the spin pseudo-gap temperature [ p^ . 
Here, we study the square lattice tight-binding model 
whose dispersion is given by e^ = 2t{coskx + cos ky) + 
4:t' cos kx cos ky + 2t"{cos2kx + cos2ky), where t, t' and 
t" denote the nearest-, next-nearest- and third- nearest- 
neighbor hopping integrals. The parameters are chosen 
as it,t',t") = (-1.0,0.17,-0.2) for YBCO and NCCO 
§, and (-1.0,0.15, -0.05) for LSCO [[|, respectively. 

First, we discuss the linear response theory for the 
TEP. In 1964, Luttinger derived a general expression for 
the TEP based on the linear response theory iQ. Al- 
though his original formula was given by a complicated 
three-particle correlation function, Jonson and Mahan 
later obtained a much simplified formula for a system 
with impurities and adiabatic phonons Pq]. In the same 
way, the formula for the TEP of a Hubbard model with 
on-site Coulomb interaction U is derived as follows Q: 

b = ^J^xx/sl (Jxx-, (1) 

L^^= $(tj + iO)/iw|^^o, (2) 

where ^(lo + iO) is given by the analytic continuation of 
the following function: 

$(l^a) = / dTer^^^u;^{TrUO)MT)), (3) 



where jk is the current operator, Tt- is a r ordering op- 
erator, and u!\ = 27riTA. Here, Oxx is the conductivity 
and — e (e > 0) is the charge of an electron. Recently, 
we performed the analytic continuation exactly as for the 
most singular terms with respect to the inverse of the 
damping rate of the quasiparticles |lq] . The result is 



XVy^x{(-l,)J\^x{f^l), 
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Jkx(e) = Wkx(e) +y^ 

k' 

J ^%iike\k'e')\G^,{e')f J^,xie'), (5) 

where Vkx{e) = (a/9A:,)(e^ + ReEk(e)), /(e) - (1 + 
exp((e— /x)/T))~^, and Gk(e) is the Green function. Here, 
Jkxi^) is the total current with the VC from the irre- 
ducible four-point vertex 7^2 j which was introduced by 
Eliashberg in ref. ||l^ for the first time. The derivation 
of cq. (0) beased on the Fermi liquid theory will be given 
elsewhere Q. 

In calculating a transport coefficient, it is known that 
the conservation laws should be satisfied to avoid un- 
physical results 111]. For this purpose, the Ward iden- 
tity between 7^2 i'^ 6q. (||) and the self-energy Sk(e) has 
to be taken into account 11^,01. Note that 722 within 
the FLEX approximation is given in ref. M. We stress 
that the relaxation time approximation cannot reproduce 
the seemingly non-Fermi liquid behaviors of Rn and the 
magnetoresistance (Ap/p): These long-standing myster- 
ies of the magnetotransport phenomena were recently 
solved by taking account of the VC for Jk in refs. ||] 
and |10|for _Rh and in refs. |13| and U^ for Ap/p, from 
the standpoint of the nearly AF Fermi liquid. Moreover, 
anomalous temperature and pressure dependences of Rn 
in K-BEDT-TTF organic superconductors are well repro- 
duced according to the same mechanism [^. Thus, it is 
important to determine whether the spin-fluctuation the- 
ory with satisfying the conservation laws is also successful 
for the thermotransport phenomena. 

As for the conductivity, we calculate ctxx based on the 
conserving approximation as explained in ref. P]. Ac- 
cording to the FLEX approximation, axx oc l/7coid be- 
cause the VC gives only a qualitative correction for axx- 
As a result, the famous T-linear resistance in high-Tc 
cuprates is reproduced for a wider range of temperatures 

Before calculating the TEP, we consider the behavior 
of eq. (0) by using the quasiparticle representation of the 
Green function, Gk(a; -t- i(5) = Zi^/{uj — ej^ -I- i7i^). Here, 
4 = ^k(e^ + ReEk(ei;)-M), 7k = - zwlmj:^{el + iO) > 0, 
and Zk — (1 — 9Sk(e)/9e)^^ is the rcnormalization factor. 
Then, eq. (H) is simplified as 



^xx — 
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where fcii [fcj^] is the momentum along [perpendicular 
to] the Fermi surface. At sufficiently low temperatures, 
eq. (||) becomes 
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Thus, Lxx given by eq. (m) is enhanced by z^. . 

When the temperature dependences of 7k(ek)j i'ka;(ek) 
and Jka;(ek) ^'^^ negligible like in a conventional Fermi liq- 
uid, then Lxx oc T"^ ll and S oa T are obtained. However, 
the temperature dependences of these functions are usu- 
ally large at higher temperatures in strongly correlated 
systems. For example, in heavy Fermion systems, a huge 
e dependence of 7(e) due to the Kondo resonance causes 
a prominent non-Fermi liquid behavior on S around the 
Kondo temperature |H,0. In the present study, we find 
that the temperature dependences of the anisotropy of 
7k (ck) ™ ^'i^- ® ^^' (0) ^^'^ primarily responsible for 
highly enhanced TEP in (under-doped) high-Tc cuprates. 
Thus, we analyze the behavior of 7k(ek) below. 
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FIG. 1. Fermi surface in YBCO/LSCO. AFBZ stands for 
the AF Brillouin zone boundary. The cold spot, at which 
ImEk(O) takes the smallest value on the Fermi surface, is lo- 
cated near (7r/2, 7r/2) [(0,7r)] in YBCO [NCCO]. 

Figure y shows the Fermi surface without interaction, 
together with the contour given by ej — ±6e (Se > 0). 
According to the ARPES measurement, the position of 
the cold spot, where 7k takes the smallest value, is near 
(7r/2, 7r/2) in YBCO and near (tt, 0) in NCCO. This "cold 
spot alternation" was first predicted using the FLEX ap- 
proximation in ref. H prior to the ARPES measurement 
P,|l. According to the analysis in ref. 0, this finding 
leads to the opposite sign of the Hall coefficient in hole- 



and electron-doped compound, although they have sim- 
ilar hole-like Fermi surfaces. Later in this article, this 
finding also explains the opposite sign of the TEP in hole- 
and electron-doped compound. 
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FIG. 2. Obtained relation between 7k = ImEk(ek) ^^id 
Ek at T = 0.02. We see that d'yk{el)/dk± > (< 0) around 
the cold spot for YBCO (NCCO), where k± is normal to the 
Fermi surface. 

Next, we examine the (k, e)-dependence of the self- 
energy Sk(e) given by the FLEX approximation. Here 
we use C/ = 8 for YBCO and C/ = 6 for NCCO. Figure | 
shows the obtained Jki^u) ^® ^k ii^ar the Fermi surface, 

e^ is the 



(tt, tt). Here 
We stress that 7k(ek) i^ 



along (0,0) -^ (7r,7r) and (0, tt) 
solution of Re{l/Gk(e;;)} = 0. 
highly asymmetric, and djit:{eD/dk± > [< 0] at the 
cold spot in YBCO [NCCO], where k± is the momentum 
normal to the Fermi surface. (Note that 7k(ek) ^ {^kl" 
and n ^ 2 if the self-energy is fc-independent.) This find- 
ing naturally explains why the sign of 5 in a hole-doped 
compound and that in an electron-doped compound are 
different, according to eq. (0) for L^x- 

Here, we analyze the origin of the asymmetric behav- 
ior of 7k (ck) i^ terms of the nearly AF Fermi liquid: In 
high-Tc cuprates, the spin propagator is well expressed 
in the following functional form: 



Xciico) = Xq/(1 + ?'(q - Q)' - iio/uJsf), 



(8) 



where Q — (tt, tt), f is the AF correlation length, and 
XQ oc u;~f^ ex ^^. According to the standard AF fluc- 
tuation theory, ^^ oc T^^ in two-dimensional systems 
[^,|l|. If we assume Wsf ^ T for simplicity, then 
7k (e) = Iml]k(e — iO) around the cold spot is approxi- 
mately given by 



7k(e) oc 



1 e{T^ + {e/T^f) 

|vk| [ l+C2(2A/fc)2]3/2' 



(9) 



within the one-loop approximation p2]. Here, Afc is 
the distance between k and the antiferromagnetic Bril- 
louin zone (AFBZ) boundary, as shown in Fig. fi^. Equa- 
tion (^ gives an analytical explanation for the numerical 



results shown in Fig. ^. (According to the FLEX ap- 
proximation, £^ < 1/Afc for n < 0.9 |9|].) By considering 
Fig. 0, eq. ^ directly indicate that (i) the inside re- 
gion of the AFBZ provides the positive contribution for 
S because d^]s_{el?)/dk±^ > in the presence of strong 
AF fluctuations. At the same time, (ii) the outside of 
the AFBZ provides the negative contribution because 
^l^.{el)/^k^ < 0. 
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FIG. 3. Temperature dependences of the TEP calculated 
with full VC's for Jk. Here, T = 0.1 corresponds to ~ 500 K. 
We put (7 = 8 for YBCO, [/ = 6 for NCCO and U = 5.5 
for LSCO. For NCCO, the VC's play a qualitatively essential 
role, and \S{n = 1.10)| > \S{n = 1.15)| is realized only when 
the VC's are taken into account. In contrast, the VC's are 
less important for YBCO and LSCO. 

In Fig. pi we show numerical results of the TEP derived 
from eq. H) according to the conserving approximation. 
The TEP without the VC for L^^, which is given by 
replacing Jka;(e) with Uk£c(e) in eq. (Q), is also plotted in 
Fig. H Here, T = 0.02 corresponds to ~ 100 K. In YBCO 



([/ = 8), dS/dT < forn = 0.80 ~ 0.92, nd S becomes 
positive at room temperatures for optimally- and under- 
doped case {n > 0.85). This result is consistent with that 
of experiments Q . We stress that S ~ a-T and a < for 
a smaller interaction {U — 3.5) as shown in Fig. H, which 
indicats the importance of the correlation effect on S. 

On the other hand, S is always negative and \S\ in- 
creases as n approaches 1 in NCCO, and dS/dT > in 
under-doped systems for higher temperatures. These re- 
sults are consistent with those of experiments |2|g] . We 
stress that the VC for Jk enhances the anomalous tem- 
perature dependence of S strongly in the case of NCCO, 
as shown in Fig. ||. Finally, dS/dT < is also repro- 
duced in the case of LSCO, and [51 for LSCO is larger 
than \S\ for YBCO or NCCO, if we compare the same 
filling cases. However, we see that the qualitative behav- 
ior of S is the same if the Fermi surface is hole-like. We 
note that S given by eq. (Q) becomes zero at T = if 
the ground state is metallic. This finding indicates that 
dS/dT should change to positive below T — 0.02 for both 
YBCO and LSCO. 
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FIG. 4. Temperature dependences of 5*1,1 and Sout, and 
R^ and 7?S"* for YBCO (n = 0.90). Here 5totai = S and 
^totai ^ ^^^ respectively. 

Finally, we discuss the following quantity to under- 
stand qualitatively the temperature and doping depen- 
dence of S for YBCO: 



Tin _ 

^xx — 



\kx\ + \ky\<iT 

E 

k 



de 



5/ 
de 



evkx (e) 



(10) 



where [ • ■ • ] is the same as that of eq. (Q). In Fig. ^ we 
show Sin = —eL™^/Taxx together with S'out = S ~ Sin- 
Apparently, S'inrout] represents the contribution from the 
inside [outside] region of the AFBZ. As we expected an- 
alytically. Sin > and S'out < is actually observed, and 
l'5'inl > I S'out I for YBCO at lower temperatures because 
the cold spots exist inside the AFBZ. In Fig. ^, we see 
a conventional behavior S'out oc T approximately. On 
the other hand. Sin ~const. is realized because l/7k(ek) 
around the cold spot becomes much asymmetric as the 
temperature decreases, as shown in Fig. ||. According to 



eq. (||), we see that d^k/dk± oc ^^Afc • 7k for ^Afc < 1, 
and 57k/9fcj_ oc (Afc)"^ • 7k for ^Afc > 1. Thus, the k- 
dependence of l/7k(ek) near the cold spot is prominent 
in the case of f ~ (Afccoid)"^, which is realized in YBCO 
for n = 0.9 at low temperatures according to our previ- 
ous study 0. Moreover, S is enhanced by z^^, which 
increases slowly as the temperature decreases. As a re- 
sult, the approximate behavior S(= Sin + Sout) oc T -I- a 
(a > 0) as well as the change of the sign of S in YBCO 
is realized. 

This situation is very contrastive to that for the Hall 
coefficient. We consider R^ which comes from the inside 
region of the AFBZ, and i?^"' = Rr - R^- As shown in 
Fig. |, \R^\ > |i?°^*| for 0.1 > T > 0.02, which means 
that i?H is almost determined by the electronic property 
at the cold spot. As a result, the simple scaling relations 
^H~ ^^ and A/9 ■ p ^ Sf' are realized in high-Tc cuprates 
ill- 

In summary, we analyzed the anomalous behavior of 
the TEP in high-Tc cuprates in a conserving manner. In 
our numerical calculation based on the spin-fluctuation 
theory, the main features of the TEP are reproduced suc- 
cessfully, at least for |1 — n| > 0.1 above the spin pseudo- 
gap temperature. The main origin is that the quasiparti- 
cle damping rate, 7k(ek)' becomes more anisotropic near 
the Fermi surface as the temperature is decreased, re- 
flecting the growth of AF fluctuations. In conclusion, 
the seemingly non-Fermi liquid behavior of the TEP is 
well understood in terms of the nearly AF Fermi liquid 
picture. In particular, the difference in the sign of the 
TEP in electron- and hole-doped systems is naturally ex- 
plained by the "cold spot alternation" mechanism. Con- 
versely, the success of the present study means that the 
AF fluctuation theory provides a reliable description for 
the excited states of quasiparticles, which determine the 
thermotransport phenomena. 

Before concluding the study, we comment that com- 
pounds with large figure of merit, Z = UxxS'^ / k {k being 
the thermal conductivity), attract great attention nowa- 
days because of their applicability in electricity genera- 
tors or refrigerators |15|. Thus, theoretical study of the 
TEP phenomena in strongly correlated systems will be- 
come much more important in the near future. 

The author is grateful to T. Saso for stimulating dis- 
cussions and comments. He also thanks K. Yamada for 
useful comments. 
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